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We investigate the possibility to detect WIMPs at the ILC in a model-independent way 
using events with single photons. The study is done with a full detector simulation of 
the LDC detector and MarlinReco. It turns out that WIMPs are observable this way at 
the ILC if their coupling to electrons is not too small. Beam polarisation can increase 
the accessible phase space significantly. 

1 Introduction 

Weakly interacting massive particles (WIMPs) are currently favoured as candidates for Dark 
Matter, which makes up about 20% percent of the total matter-energy content of the uni- 
verse. Although the lightest supersymmetric particle, which is the lightest neutralino in 
many SUSY models, would make a very good WIMP if SUSY with R-Parity conservation 
is realised in nature, most extensions of the Standard Model of particle physics contain 
WIMPs. In fact the more general requirement is that some new conserved quantum number 
makes the lightest of the new particles stable. The International Linear Collider (ILC) offers 
the possibility to look for WIMPs in a rather model-independent way, which first has been 
pointed out in [2]. Here, the expected sensitivity to such a WIMP signal, the achievable 
mass resolution and the influence of beam polarisation are studied using a full detector 
simulation. 

Assuming that the cosmic relic density of WIMPs is determined by pair annihilation of 
WIMPs, and that an unknown branching fraction Kg of these annihilations proceeds into an 
e"'"e~ paiv^XX — > e~^e~), crossing relations can be used to derive an expected cross-section 
for the reverse process, i.e. e~^e~ XX . This cross-section contains as free parameters: 

• the e+e^ branching fraction Ke 

• the mass of the WIMP Mx 

• the spin of the WIMP Sx 

• the angular momentum of the annihilation's dominant partial wave J. 

In order that this process be observable in a collider detector, where the WIMPs them- 
selves leave no signature, an additional photon from initial state radiation is required: 
e+e" XX-f. 

The main Standard Model background process for this reaction is neutrino pair pro- 
duction, again with an ISR photon: e+e^ —* vvj. At energies significantly above the Z'^ 
pole, this reaction is dominantly mediated by i-channel VF-exchange and can thus be re- 
duced significantly by choosing the appropriate polarisations for the electron and position 
beams. Therefore we consider three possible scenarios for the helicity structure of the WIMP 
coupling to electrons: 
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• the same as the SM charged current weak interaction, i.e. only K(e^e^) is nonzero 

• parity and hehcity conserving, i.e. ^(e^ej^) = ^(ejje^) 

• opposite to SM charged current weak interaction, i.e. only = ^{^r^^,) is nonzero. 

Especially in the last case a significant enhancement of the signal over background ratio 
is expected. 

2 Software and Reconstrucion Tools 

The neutrino background has been generated with NUNUGPV [3], which is a dedicated 
generator for neutrino pair production with up to three photons e~'"e~ — > 1.2 • 10^ 

events with at least one photon with an energy 8 GcV < < 250 GeV and a polar angle 
9.y 15° < 9^ < 165° were generated at a center-of-mass energy of ^/s = 500, corresponding to 
an integrated luminosity of 500 fb~^. These events serve not only to describe the irreducible 
Standard Model background, but arc also reweighted w.r.t. energy and polar angle of the 
photon according to the WIMP cross-section. The benefit of the is method is that in this 
way a signal sample can be obtained for all mass and spin hypotheses to be tested without 
applying the detector simulation again, which reduces the processing time tremendously. 

This sample has been subjected to the full LDC detector simulation, using the detector 
model LDCOlSc with a 4 T magnetic field and Mokka 6.1 [4]. The events were then recon- 
structed with MarlinReco [5], using the WOLF [6j algorithm for particle flow and a simple 
selection demanding E^ < 10 GeV and 20° < 6^ < 180°. For the energy resolution studies 
described in the following, also an angular match to the generated photon was required. 

3 Energy Resolution Studies 

The infiucncc of the detector's energy resolution has been studied at two levels: besides 
the full reconstruction, the so called cheated reconstruction makes use of the simulation 
information for associating calorimeter clusters with particles. This way the pure detector 
resolution and calibration can be disentangled from confusion effects, noise and so on. Fig- 
ure[l]shows the generated (left) as well as the cheated and fully reconstructed (right) energy 
spectra of the most energetic photon candidate of the events. 

While the peak at E^ = 240 GeV from radiative returns to the Z° resonance is extremely 
clear at the generator level, it his heavily smeared after reconstruction, even for the cheated 
case. The full reconstruction yields even less photons at high energies. The reason for this 
effect can be seen in the left plot of figure [21 which shows the mean number of photon 
candiates per generated photon as a function of the generated photon's energy. 

It is clearly visible that at high photon energies the WOLF algorithm tends to split 
clusters stemming from one photon into several photon candidates. Therefore we apply 
a merging procedure to neighboring photon candidates, after which the fully reconstructed 
photon energy distribution is practicaly identical to the cheated spectrum shown in the right 
hand plot of figure [TJ 

The energy resolution obtained after the recombination procedure is shown in right plot 
of figure [2] as a function of 1/^/E. With being roughly constant at about 6%, it seems to 
be significantly worse than the 14.4%\/E + 0.5% aimed for at the ILC. This effect has been 
meanwhile tracked down to an imperfect calibration, which means that the results presented 
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Figure 1: Energy spectra of the event's most energetic photon at generator level (left) as 
well as cheated and fully reconstructed (right) 



in the next section are expected to improve further once a better calibration procedure is 
applied. 
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Figure 2: Left: Mean number of reconstructed photon candidates per generated photon vs 
the generated photon's energy. 

Right: Photon energy resolution after the recombination procedure vs . 



4 Preliminary Analysis Results 

Currently, the following scenarios have been investigated: 

• WIMP spin: P-wave annihilation (J=l) for Sx — 1 and Sx ~ \ 

• WIMP couphngs: K(e£eJ) > 0, K(e]^eJ) > and ^(e^ej) = K(e^eJ) > 

• polarisation: unpolarised beams, e~ polarisation only {P^- ~ 0.8) and additional e 
polarisation (Pg- = 0.8 and Pe+ = 0.6). 
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Observation reach 

For each combination of these parameters, the reach of the ILC with an integrated himinosity 
of 500 fb~^ at = 500 GeV for a 3cr observation of WIMPs has been determined as a 
function of the WIMP mass. Due to the high irreducible background from Standard Model 
neutrino production, the sensitivity has been obtained statistically by using fractional event 
counting [7] as implemented in the RooT class TLimit. 
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Figure 3: 3cr observation reach of the ILC for a Spin-1 WIMP in terms of WIMP mass and 
for three different assumptions on the chirality of the electron- WIMP coupling, see text. Full 
line: P^- = = 0, dotted line: P,- = 0.8, P<.+ = 0, dashed line : P^- = 0.8, P^+ = 0.6. 
Regions above the curves are accessible. 

Figure [3] shows the expected ILC sensitivity for Spin-1 WIMPs in terms of the mini- 
mal observable branching fraction to electrons Kg as a function of the WIMP mass. The 
leftmost plot shows the case where the WIMPs couple only to lefthanded electrons and 
righthanded positrons (^(e^ej)), the middle plot shows the parity and helicity conserving 
case (^(e^ej) = K(ej^eJ), while the right plot is dedicated to the case that the WIMPs 
couple to righthanded electrons and lefthanded positrons ('t(e]^e^). The regions above the 
curves are accessible, where the full line gives the result for unpolarised beams, the dotted 
line for Pg- = 0.8 and the dashed line for Pg- = 0.8 and Pg+ = 0.6. In the latter two 
coupling scenarios polarised beams increase the reach significantly, especially the additional 
positron polarisation increases the accessible range in Kg by about a factor of 2. Figure |4] 
shows the same for a Spin-^ WIMP. Here the sensitivity is somewhat worse, but again beam 
polarisation extends the observable part of the parameter space significantly. 




Figure 4: 3tT observation reach of the ILC for a Spin-i WIMP in terms of WIMP mass and Kg 
for three different assumptions on the chirality of the electron- WIMP coupling, see text. Full 
line: Fg- = Pg+ = 0, dotted line: Pg- = 0.8, Pg+ = 0, dashed line : Pg- = 0.8, Pg+ = 0.6. 
Regions above the curves are accessible. 
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Mass resolution 



polarised beams 



If WIMPs are observed at the ILC, their mass can be determined from the recoil mass 
distribution of the photons: 

M2 = s - 2^E^ (1) 

Figure [H] shows an example for the recoil mass distribution for a 150 GcV Spin-1 WIMP 
with both beams polarised. The WIMP signal shown in dark grey kicks in at Afiocoii = 
316 GeV. 

From this distribution, the WIMP mass 
can be reconstructed for example by a tem- 
plate method. For this procedure, only 
200 fb~^ of the available MC sample have 
been analysed as dataset, the rest is used for 
the templates. Figure [5] shows the obtained 
Ax^ as function of the reconstructed WIMP 
mass for a 150 GeV Spin-1 WIMP for = 
0.3. The left plot shows the helicity and par- 
ity conserving case, the right plot the case 
that the WIMPs couple to righthanded elec- 
trons and lefthanded positrons (^(e^^ej). 
Again the full line gives the result for unpo- 
larised beams, the dotted line for Pg- = 0.8 
and the dashed line for P^- ~ 0.8 and 
Pe+ = 0.6. Without any beam polarisation, 

the mass resolution is about 4 GeV, which is reduced to about 1.2 GeV by switching on the 
electron polarisation. Positron polarisation improves the resolution by another factor 2 to 
about 0.6 GeV. 




100 150 200 290 300 390 400 490 900 
MrecollEGeV] 



Figure 5: Recoil mass 
150 GeV Spin-1 WIMP 



distribution for a 
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Figure 6: Ax^ from mass determination by a template method for a Spin-1 WIMP with 
Mx=150 GeV. Left: parity and helicity conserving couplings, right:K(e]^eJ). Full line: 
^ ^ 0, dotted line: P^- = 0.8, Pe+ = 0, dashed line : P^- = 0.8, P,+ = 0.6. 

Figure [7] shows the analoguos results for a 180 GeV Spin-i WIMP. As for the observation 
reach, the situation is slightly worse than in the Spin-1 case, but again the use of beam 
polarisation leads to a significant gain in resolution. 



LCWS/ILC2007 




[GeV] [GeV] 

Figure 7: Ax^ from mass determination by a template method for a Spin--^ WIMP with 
Mx=150 GeV. Left: parity and heUcity conserving couphngs, right:«;(e^ej). FuU hne: 
= = 0, dotted hne: P^- = 0.8, Pe+ = 0, dashed Hne : P^- = 0.8, Pe+ = 0.6. 

5 Conclusions and Outlook 

The study of model-independent WIMP production at the ILC presented here is one of the 
first examples of analyses performed using the full simulation of the LDC detector. A recon- 
struction for these fully simulated events exists and usable for analysis, but needs further 
improvements in parallel to the optimization of the detector concept. Already with the cur- 
rent level of sophistication this study shows that there is a good chance to detect WIMPs in 
this model-independent way at the ILC and to measure their mass with a precision of about 
1 GeV. Both the range in phase space as well as the mass resolution improve significantly 
when polarised beams are assumed. Typically the use of 80% electron polarisation gives 
improvements of a factor of two over unpolarised beams, whereas an additional positron 
polarisation of 60% yields another factor of two. 

The results presented here will be improved in the near future by applying a more 
appropriate detector calibration and by using other particle flow algorithms and photon 
finders. Furthermore reducible backgrounds as well as beamstrahlung have to be included 
in the study, before finally different detector concepts can be compared. 
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